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Abstract 

Recent results on sub jet distributions from ZEUS are presented. The measured nor¬ 
malised cross sections were used to study the pattern of parton radiation. The compari¬ 
son of the measurements with leading-logarithm parton-shower Monte Carlo models and 
perturbative QCD calculations shows a good agreement between data and predictions. 
Results on event-shape mean values and distributions are also presented. These measure¬ 
ments were used to test the predictions of the power-correction model for hadronisation. 
A universal value, within 10%, of the effective parameter do of the model was obtained. 


1 Introduction 

At lowest-order QCD, the diagrams that contribute to neutral current (NC) deep inelastic ep 
scattering (DIS) at HERA are the boson-gluon fusion (BGF) (Vg —■► qq, where V = 7 * or 
Z°, Fig. Eh) and QCD-Compton (QCDC) (Vq —> qg, Fig. Eh) processes. Figure Et shows an 
example of a higher-order diagram. Any observable can be expressed as the convolution of the 
parton densities in the proton, f a , times the matrix elements, 

^2 fa(x, Rf) ® Matrix Elements, 

a 

where r^ is the factorisation scale. The matrix elements, in addition to x and (jlf, depend 
also on the strong coupling constant, a s , and the renormalisation scale Rr. The sum runs over 
all type of partons a. Thus, in the regions of phase space where the parton densities are well 
constrained, measurements of e.g. jet cross sections can be used to perform tests of perturbative 
QCD (pQCD) and determinations of a s . 

The hadronic final state in NC DIS can also be used to study the pattern of parton radiation 
by means e.g. of subjets inside jets. Subjets observables are calculable in pQCD and so 
they provide stringent tests of the theory. On the other hand, the hadronisation process, a 
non-perturbative effect, can be studied by means of the event shapes. Recent developments 
on the model of power-law corrections |Tj have prompted revived interest in understanding 
hadronisation within the framework of pQCD. 

*Talk given at HEP2005 International Europhysics Conference on High Energy Physics, Lisbon, Portugal, 
July 21 st - 27 th , 2005. 
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Figure 1: Examples of Feynman diagrams for deep inelastic ep scattering processes at lowest- 
order QCD: (a) boson-gluon fusion and (c) QCD Compton, (c) Example of a higher-order 
diagram. 


2 Subjet distributions 

The investigation of the internal structure of jets gives insight into the transition between 
a parton produced in a hard process and the experimentally observable jet of hadrons. At 
sufficiently high jet transverse energy, , where the effects of fragmentation can be neglected, 
the jet structure can be calculated perturbatively. The lowest non-trivial-order contribution to 
the jet substructure is given by 0(a s ) calculations for NC DIS in the laboratory (LAB) frame. 
Next-to-leading-order (NLO) calculations of jet substructure can be obtained in the LAB frame 
since, in such a case, it is possible to have three partons inside one jet. 

The kr cluster algorithm |2] was used in the longitudinally invariant inclusive mode t 3j to define 
jets in the hadronic final state. The internal structure of the jets can be studied by means of 
the subjet topology. Subjets were resolved within a jet by considering all particles associated 
with the jet and repeating the application of the k t algorithm until, for every pair of particles 
i and j the quantity 


dij = min (E Tti , E T>j ) 2 • [(^ - rjj) 2 + (<& - fa) 2 ], 

where E T rji and 0* are the transverse energy, pseudorapidity and azimuth of particle i, 
respectively, was greater than d cut = y cnt • (E^) 2 . All remaining clusters were called subjets. 
The subjet multiplicity depends upon the value chosen for the resolution parameter y cut . 

The pattern of QCD radiation from a primary parton has been studied [lj by measuring nor¬ 
malised cross sections as a function of subjet observables: the ratio between the subjet trans¬ 
verse energy and that of the jet, E^/Ej?, the difference between the subjet pseudorapidity 
(azimuth) and that of the jet, r] sh i — rf et (|0 sb l —0i et |), and cc sb j, the angle, as viewed from the jet 
centre, between the highest transverse energy subjet and the beam line in the pseudorapidity- 
azimuth plane. The measurements were done for Q 2 > 125 GeV 2 , where Q 2 is the momentum 
transfer. Jets of E^ > 14 GeV and —1 < rf et < 2.5 were selected. The final sample consisted 
of those jets which had two subjets for y cut = 0.05. 

The 0(a s ) and 0(a 2 ) QCD calculations used to compare with the data are based on the 
program DlSENT |3]. For these calculations, the number of flavours was set to five; the renor¬ 
malisation and factorisation scales were both set to Hr — /if — Q] was calculated at two 
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loops using = 220 MeV, which corresponds to a s (M z ) = 0.1175. The MRST99 [Hj pararn- 
eterisations of the proton parton density functions (PDFs) were used. 

The cross-section (1 /a)(da/d(E^/E^)) is presented in Fig. The distribution of the fraction 
of transverse energy contains two entries per jet and is symmetric with respect to = 0.5 

by construction. The data distribution has a peak at E^/E^ = 0.5, which shows that the 
two subjets tend to have similar transverse energies. The distribution for the difference in 
pseudorapidity is shown in Fig. EJo and also has two entries per jet. The measured cross section 
has a two-peak asymmetric structure, with a dip at rj sbi — ?/l et ~ 0, which shows that the subjets 
cannot be reconstructed too close together. FigureEh presents the normalised cross section as a 
function of 0 st> j — <^ et |. There are two entries per jet in this distribution. The data distribution 
has a peak at |0 sb j — <^ et j = 0.2 — 0.3; the suppression at |</> sb -i — <^ et | ~ 0 comes also from the 
fact that the subjets cannot be resolved when they are too close together. The distribution as 
a function of a sb -j (one entry per jet) increases as a sb -i increases (see Fig. Et>)- This shows that 
the highest transverse energy subjet tends to be in the rear direction. This is consistent with 
the asymmetric peaks observed in the rj sh i — r]^ et distribution. 
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Figure 2: Normalised differential subjet cross sections as functions of (a) E^/E^ and (b) 
^sbj _ ^jet. p Qr comparison, the predictions of Disent at LO (dashed histograms) and NLO 
(solid histograms) are included. 

The fixed-order QCD calculations are compared to the data in Figs. El and El The QCD 
predictions give a good description of the data in shape, within 10%. This shows that the 
mechanism driving the subjet topology are the q —> qg and g —> qq subprocesses as implemented 
in the pQCD calculations. 

To study in more detail the pattern of parton radiation, the predictions of quark- and gluon- 
induced processes are compared separately with the data in Fig. El The NLO calculations 
predict that the two-subjet rate is dominated by quark-induced processes: the relative contri- 
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Figure 3: Normalised differential subjet cross sections as functions of (a) |0 st| j — <^ et | and (b) 
o sb j . For comparison, the predictions of Disent at LO (dashed histograms) and NLO (solid 
histograms) are included. 


bution of quark- (gluon-) induced processes is 82% (18%). The predictions for these two types 
of processes are different: in quark-induced processes, the two subjets have more similar trans¬ 
verse energies (Fig. 0b) and are closer to each other (Figs. 0Jo and 0b) than in gluon-induced 
processes. The comparison with the measurements shows that the data are better described 
by the calculations for jets arising from a qg pair than those coming from a qq pair. 


3 Event shapes 

Event-shape observables are particularly sensitive to the details of the non-perturbative effects 
of hadronisation and can be used to test the models for these effects. In this type of analysis, 
the data are compared to model predictions which combine NLO calculations and the theoret¬ 
ical expectations of the power-corrections (PC) model, which is characterised by an effective 
coupling «o- The total prediction for any event-shape observable is then given by the sum of the 
perturbative and PC predictions. Previous results supported the concept of power corrections 
in the approach of Dokshitzer et ah, but a large spread of the results suggested that higher- 
order corrections were needed. Now, resummed next-to-leading-logarithm (NLL) calculations 
matched to NLO are available and so it is possible to study event-shape mean values as well as 
distributions. 

The event-shape observables studied are thrust, T, broadening, B, the C parameter and the 
jet mass, M, 
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Figure 4: Normalised differential subjet cross sections as functions of (a) E b ^/E ^ Xl , (b) r) sh i—rji et 
and (c) |0 sb -i — 0J et |. For comparison, the NLO predictions for quark- (solid histograms) and 
gluon-splitting (dashed histograms) are included. 

2 (E,:g.) 2 -IZ,p.l 2 
2(Eilftl) 2 ' PET,) 2 ’ 

where n specifies the thrust or 7 axis, 6 t] is the angle between two final-state particles i and j, 
and Pi ( Ei ) is the momentum (energy) of particle i. 

A suitable frame in which to study event shapes at HERA is the Breit frame since in this frame, 
the separation between the current jet and the proton remnant is maximal. The event-shape 
variables are reconstructed for all the particles in the current hemisphere of the Breit frame. 

Measurements of event-shape means have been made (I) as a function of Q for each observable 
in the kinematic region given by 80 < Q 2 < 2 • 10 4 GeV 2 and 0.0024 < x < 0.6 (see Fig. EH). 
Predictions consisting of NLO + PC calculations have been fitted to the data, leaving a s 
and do as free parameters. Each observable was fitted separately. The NLO predictions were 
calculated using the program Disaster++ 5j. The proton PDFs have been parameterised 
using the CTEQ4M |2| sets. A reasonable fit is obtained for all the event-shape observables 
within the Q 2 range studied. 

Figure EH shows the event-shape differential distributions [7j for some of the observables in 
different regions of Q. In this case, the fit was done using NLO + NLL + PC. The NLO 
predictions have been calculated using the Disaster++ and Dispatch PI programs with 
the MRST99 proton PDF sets. The PC and matched NLL predictions were calculated using 
the Disresum ma package. In Disresum, the power correction is applied as a shift of the 
distribution, which has the same functional form as the power correction for the mean. For B 1 , 
there is in addition a change in shape. The range of the fit for each observable was restricted to 
the regions where the predictions were valid. A reasonable fit is obtained for all the event-shape 
observables within the restricted ranges studied. 

The extracted values of a s and do from the means and differential distributions are shown in 
Figs. Ek and Et>, respectively. It is possible to obtain a universal value for cto of 0.45 at the 
10% level, except for T 7 (means) and C parameter (distributions). The extracted values of do 
and a s (Mz) from the means show a dispersion that could be due to higher-order terms. The 
extracted values of a s (Mz) from the distributions are consistent with the world average. 
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Figure 5: The means of event-shape observables fitted with NLO + PC. 

Measurements of an event-shape observable sensitive to higher-order effects have been made [7j. 
The out-of-plane momentum, /Pout = |P? ut |> which is defined as the energy flow out of the 

plane defined by the proton direction and the axis which maximises the thrust, is sensitive 
both to perturbative and non-perturbative contributions. The lowest non-trivial contribution 
comes from non-perturbative effects or from 0(a 2 ) contributions. Figured shows the measure¬ 
ments of A out /Q in two ranges of Q 2 . The data are well described by the predictions of the 
leading-logarithm parton-shower model of Lepto HU. This constitutes the first comparison 
of calculations which include leading-order + NLL + PC with the data, in the liigh-Q 2 range 
only. The description of the data by the prediction is reasonable, but a more precise test of the 
model needs higher-order calculations. 


4 Summary 

Subjet normalised cross sections have been measured in NC DIS using 81.7 pb -1 of data col¬ 
lected with the ZEUS detector at HERA with a centre-of-mass energy of 318 GeV. A reasonable 
description of the data is obtained by the QCD predictions. This means that the pattern of 
parton radiation as implemented in a NLO calculation reproduces the behaviour of the data. 
In addition, the data are well described by the calculations for jets arising from the splitting of 
a quark into a quark-gluon pair. 

Event-shape means and distributions have been measured in NC DIS using 82.2 pb -1 of data 
collected with the ZEUS detector at HERA with a centre-of-mass energy of 318 GeV. Calcula¬ 
tions including NLO + PC, and including resummed NLL predictions matched to NLO for the 
differential distributions, give a reasonable description of the data. The extracted values of the 
power-correction parameter, a 0 ) are consistent within 10%. The extracted values of the strong 
coupling constant are consistent with the world average. However, more theoretical input is 
needed to fully exploit the potential of these measurements. 
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Figure 6: Event-shape distributions as functions of (a) B 7 , (b) M 2 and (c) C, fitted with NLL 
resummed calculations matched to NLO + PC. 
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Figure 7: Extracted a s (M z ) and a 0 parameters from (a) mean values and (b) distributions, (c) 
Distribution of A'qut /Q compared with Lepto and leading-order + NLL + PC calculation. 
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